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Abstract Glued timber products have an extensive range 
of applications in construction. In this work a Microfocus 
X-ray Computed Tomography method was developed to in¬ 
spect gluing defects in timber samples and was applied suc¬ 
cessfully on experimental data. The bonding plane was seg¬ 
mented into glued and non-glued regions and imaged with 
5 mm resolution. Moreover, the gap topology between tim¬ 
ber lamellas was precisely characterised. A limited-angle re¬ 
construction with anisotropic frame binning together with a 
specific glue line readout method efficiently filters out un¬ 
desired wood structure highlighting the information of the 
adhesive joint. This method imposes limitations on the size 
of the specimen in only one dimension. The presence and 
absence of glue could be detected for glue line thicknesses 
over 50 pm and air gaps larger than 150 pm could be char¬ 
acterised. Several information reduction approaches were 
combined in the reconstruction process to implement the as¬ 
sessment of a 100 x 100 mm 2 bonding plane in less than 

40 s. 


zur Beurteilung von Fehlverklebungen in Holzbauteilen 
entwickelt und in experimentellen Versuchen erfolgreich 
iiberpriift. Fehlverklebungen und die Topologie der Kleb- 
fugen wurden mit einer raumlichen Auflosung von 5 mm 
abgebildet. Verfahrensspezifische Bildrekonstruktions- und 
Klebfugen-Auslesealgorithmen mit eingeschranktem Auf- 
nahmewinkelbereich filtern unerwiinschte Holzstrukturen 
heraus und heben die Leimschicht hervor. Die GroBe der 
Holzbauteile ist bei diesem Messverfahren nur in einer Di¬ 
mension beschrankt. Die An- und Abwesenheit von Kleb- 
stoff lasst sich fur Leimfugen dicker als 50 pm nachwei- 
sen; Luftspalten groBer als 150 pm sind feststellbar. Die 
Informations-Reduktionsverfahren ermoglichten die Beur¬ 
teilung von 100 x 100 mm 2 verleimten Holzoberflachen in 
einer Messzeit von weniger als 40 s. 


1 Introduction 


The conventional forms of sawn wood in construction have 
been replaced during recent years in favour of highly en¬ 
gineered glued timber products. The bonding quality needs 
to be assessed in order to maximise the life cycle of the 
constructions. Current European standardised methods rely 
on visual inspection or mechanical tests on selected speci¬ 
mens (EN 391:2001; EN 14080:2005). The development of 

a non-destructive glue line assessment procedure is an open 
research question. Infrared thermography detects changes 
in thermal conductivity or thermo-mechanics in delami¬ 
nated wood veneer of a few millimetres thickness (Berglind 
and Dillenz 2003; Choi et al. 2008). Ultrasonics is an ir¬ 
replaceable tool for adhesion testing due to its high sen¬ 
sitivity to delaminated interfaces and relatively low cost 
(Kabir et al. 2002; Maeva et al. 2004; Solodov et al. 2004; 
Dill-Langer et al. 2005; Hasenstab 2006). The main draw¬ 
backs are poor, coupling pressure dependent reproducibil¬ 
ity and the difficulties of performing imaging with sufficient 


Detektion von Fehlverklebungen in Holzbauteilen 
mittels Mikrofokus Rontgen-Computertomographie 


Zusammenfassung Holzverklebungen gewinnen zuneh- 
mend an Bedeutung. Im Rahmen dieser Arbeit wurde ein 
Mikrofokus Rontgen-Computertomographie Messverfahren 
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spatial resolution. Recently an air-coupled ultrasound sys¬ 
tem has been developed that is not limited by these short¬ 
comings (Sanabria et al. 2010a). 

Computed tomography (CT) allows for a full 3D recon¬ 
struction of the interior of a sample regardless of its exterior 
geometry. In the field of wood science there exist, amongst 
others, applications for imaging, segmentation and quantifi¬ 
cation of bark, knots, heartwood/sapwood boundary, spiral 
grain, moisture transport, decay, cracks, voids and year ring 
structure (Niemz et al. 1997; Schmoldt et al. 2000; Oja 2000; 
Sepulveda et al. 2002; Bucur 2003; Rinnhofer et al. 2003; 
Espinoza et al. 2005; Bhandarkar et al. 2006; Longuetaud 
et al. 2007; Scheepers et al. 2007; Entacher et al. 2007; 
Mannes et al. 2009a). However, few works on CT glue line 
assessment in timber components are available. Commercial 
X-ray density profiler systems perform a global assessment 
of the glue line on small specimens (50 x 50 mm 2 ) with 
a resolution down to 10 (dm (QDP-01X, Quintek Measure¬ 
ment Systems, Knoxville, TN, USA; Dense-Lab X, Elec¬ 
tronic Wood Systems Int., Beaverton, OR, USA; DAX 5000, 
Fagus-GreCon Greten GmbH & Co. KG, Alfeld, Germany). 
Sirr and Waddle (1999) performed X-ray CT imaging of 
animal glue repair in Stradivarius violins; Morigi et al. 
(2007) visualised stucco fillings in paintings on wood; to- 
moslices showing the presence of glue in finger joints have 
been obtained using X-rays (Niemz et al. 1999; Hu and 
Gagnon 2007) and fast neutron CT (Mannes et al. 2009b; 
Osterloh et al. 2008). 

In this work the application of Microfocus Computed 
Tomography (|dCT) to glue line inspection of timber lami¬ 
nates is shown. A limited-angle inspection procedure is used 
which greatly reduces data acquisition and reconstruction 
times, imposes limitations to only one of the dimensions of 
the specimen and eliminates undesired wood structure from 
the tomoslices. Next, specific glue line readout and segmen¬ 
tation algorithms are presented, which retrieve a spatial map 
of the areas with and without adhesive within the bonding 
plane, providing also information on the topology of the gap 
between timber lamellas. Finally, the method is optimised 
by minimising the amount of input data required to perform 
a satisfactory glue line assessment. 


Purbond AG, Sempach Station, Switzerland), was applied 
to one side of the lamellas with an amount of 200 gm -2 . 
Only two areas of about 100 x 30 mm 2 on both left and 
right edges of the sample were glued. The lamellas were 
pressed together in a hydraulic press for 3 hours at a pressure 
of 0.8 Nmm -2 . Normalised climatic conditions (T = 20°C 
and RH = 65%) were used for storage. 

The cross-section of the samples is approximately per¬ 
pendicular to the growth direction of the stem (R-T plane), 
the ring angle varying between 90° (year rings perpendic¬ 
ular to the glue joint) and 45°. Aluminium sheets of 500 
and 1000 jam thickness were introduced in the central posi¬ 
tion of the specimens (250 mm from the edges) in order to 
achieve a defined air gap thickness distribution between non- 
glued wood lamellas (Fig. 1). The gap was measured from 
the external surfaces of the sample using an optical scan¬ 
ner with 600 dpi digital resolution. The thickness of the glue 
line and the air gap in the vicinity of the transition between 
glued and non-glued areas was characterised with a high res¬ 
olution stereomicroscope with a zoom range of 7.8 x-160x 
(M205C, Leica Microsystems GmbH, Wetzlar, Germany). 
After the investigations, the samples were broken up and the 
profile of the transition between glued and non-glued areas 
was highlighted with a pencil and recorded with the optical 
scanner. 


2.2 X-ray data acquisition 


X-ray measurements were performed with a laboratory cone 
beam jaCT system optimised for non-destructive testing 
(Fig. 1). The source is a 160 kV microfocus X-ray tube 
(TXD9160, Viscom AG, Hannover, Germany) with a trans¬ 
mission target consisting of diamond that was coated with a 
3 jam tungsten layer. The tube was operated at an accelera¬ 
tion voltage of 40 kV and a nominal power of 30 W to ob¬ 
tain a focal spot size of approximately 50 jam with sufficient 
signal-to-noise ratio. The detector is an X-ray flat panel sen¬ 
sor (7942CA-02, Hamamatsu Photonics K. K, Hamamatsu 
City, Japan) with a sensitive area of 120 x 120 mm 2 and a 
pixel size of 50 jam, consisting of a CsI:Tl scintillator plate 
with needle structure. An exposure time of 1 s was used and 
the aluminium front-face of the detector was replaced with 
black paper so as to increase the efficiency for the main part 
of the tube spectrum (around 25 keV). 

The glued timber samples were placed on a motorised 
low runout rotary stage located between the X-ray tube and 
the detector and rotated with specific angular steps of 0.15° 
or 0.2°. For each angle the logarithmic X-ray attenuation 
between the source and the pixels of the detector (w/, vj) 
was recorded as a projection or frame P$> k (,w z -, vj) of 12 bits 
grey values; light grey corresponds to high attenuation. At 
cp£ =0°, the glue line plane is perpendicular to the detector 
and aligned with the main axis of the cone beam. Fimited 


2 Materials and methods 


2.1 Sample preparation and optic measurements 


A total of six glued timber samples were taken from the 
available pallet used in the air-coupled ultrasound investi¬ 
gations (Sanabria et al. 2010a). They each consisted of two 
common spruce (Picea abies Karst.) lamellas of dimensions 
500 x 100 x 5 mm 3 glued together to form a glued tim¬ 
ber object of 500 x 100 x 10 mm 3 . The adhesive, a one- 
component moisture-curing polyurethane resin (HB 110, 
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Fig. 1 (iCT measurement setup 
and definition of sample and 
X-ray detector coordinate 
systems 

Abb. 1 pCT Messaufbau und 
Koordinatensysteme fur 
Holzproben und 
Rontgendetektor 


X-ray detector 


D sd = 620 mm 


D so = 504 mm 




X-ray source 
(cone beam) 


Glued timber sample 


] were obtained by consider- 


computer with a 2.83 GHz Intel® Xeon® E5440 proces¬ 
sor and 32.0 GB of RAM memory. A Shepp-Logan filter 
was applied to P<$ k (m*, vj) for grey level contrast optimisa¬ 
tion. The limited-angle backprojection was implemented by 
summing Ao scaled and filtered frames P® k (u , v ) using a 
bilinear interpolation procedure: 


angular ranges [—O 
ing only a limited number of equispaced projections in the 
reconstruction algorithm. The distances between source and 


O 


max> 


max 


detector (D s d = 
rotation axis (D 
largest possible field of view in the detector (112 mm) with 
a sample-related pixel size of 41 jam. Two jaCT measure¬ 
ments A and B per sample were performed, for each of them 
a volume of 100 x 100 x 10 mm 3 containing glued and non- 
glued timber was reconstructed. 


620 mm) and between source and sample 

504 mm) were chosen to achieve the 
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D sd xi sin + y m cos 
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2.3 Limited-angle reconstruction 


Dsd Zn 


v(xi,y m ,Zn, ®k) 


Figure 2 summarises the main signal processing steps ap¬ 
plied to the measured jaCT data. The theoretical interpre¬ 
tation of the limited-angle method is based on the Fourier 
Slice Theorem (Kak and Slaney 1988). For an inspec¬ 
tion limited to the angular interval [—O 
Fourier transform of the CT tomoslices is only defined for 
[tt/2 - O 

angles outside [—O 
information is restricted to a quasi-horizontal angular range 
(Fig. 3). Anatomical wood features such as the year ring 
structure are distributed over the full angular range and are 
eliminated in the reconstruction process. A modification of 
the FDK cone beam filtered backprojection reconstruction 
algorithm (Feldkamp et al. 1984; Kak and Slaney 1988) was 
implemented with single-threaded Matlab® software in a 
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where /jl (v/, y m , z n ) is the reconstructed attenuation coeffi¬ 
cient. Hereafter, the discretization subscripts /, j, k, /, m and 
n are dropped in order to simplify notation. The amount of 
reconstruction data was minimised with a qualitative anal¬ 
ysis of the glue line segmentation results. Larger detector 
pixel sizes were implemented with overlapped averaging of 
Po(w, v ) followed by anisotropic binning; the reconstruc¬ 
tion step in y equals the binned pixel size in u , the steps 
in v and z correspond to the one in v . Local tomography al¬ 
lowed for further data reduction by cropping P$(u, v ) in u 
to the effective field of view FOV e necessary to reconstruct 


T* max ] the 2D 
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]; therefore, image features with 
] are filtered out. The glue line 
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pCT projections P m {u,v) 
Occurred hypotheses H = {h a ,h g \ 


Estimation of ju(x,y,z) with limited-angle 
FDK cone beam reconstruction 


Extrema computation ju ext (x,z,h) 

Mexl ( h a) = min M 

yeDs 


Mext {fig ) 


max // 
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h a : false, h 
h a : false, h 


Greyscale 

assessmen 
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I Binary 
assessment 
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Estimation of iv 
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false 


H = h 
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ju>t 
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Fig. 2 Flow diagram of glue line assessment algorithm 
Abb. 2 Flussdiagramm des Klebfugen-Auslesealgorithmus 


a volume containing the glue line plane. FOV e was calcu¬ 
lated with taking only into account ray trajectories between 
source and detector that go through the volume: 


for wood species and adhesives used in construction (Dunky 
and Niemz 2002). In the case of non-glued material, an air 
gap with a reduced attenuation coefficient is expected. Grey 
level profiles were drawn across the glue line for each bond¬ 
ing plane position (Fig. 3). The absolute extrema of each 
profile n ex t{x, z, h) were calculated; for h g and h a absolute 
minima and maxima were searched, respectively: 

z, h a ) = min /z(v, y, z) 

yeD s 

fi ext (x,z,h g ) = ma x[i(x,y,z) 

y^D s 

where D s is an interval where the glue line is certain to be 
found, typically |Z) 5 1 = §. The mean fi w and standard devi¬ 
ation S/a of the grey level in the wood pixels around the glue 
line were also computed for each profile. A preliminary de¬ 
cision was taken on h a or h g by locally testing the following 
inequalities: 

\fX e xt(h)-fl 

Ipextiha) Mw\h < I Fextihg) 

where k and X are two empirical constants. The first inequal¬ 
ity checks the contrast of the extrema with respect to back¬ 
ground noise caused by wood heterogeneity (typically k = 2 
is sufficient). The second one determines which hypothesis 
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where w is the width of the sample in direction v at O = 0° 
and § is the thickness reconstructed across the glue line. For 
the setup with u ! = 100 mm and § = 10 mm the glue line 
assessment algorithms performed well. 


( 3 ) 


2.4 Glue line readout algorithm and modified MAP 

segmentation 


The segmentation of glued and non-glued areas was ad¬ 
dressed as a binary hypothesis test, with h g the hypothesis 
of presence of glue and h a the hypothesis of occurrence of 
an air gap between timber lamellas. First a specific glue line 
readout applied to the tomoslices (Fig. 2) reduced the 3D 
tomogram /x(v, y, z) to a summarised 2D greyscale image 
of the bonding plane jl(x,z ). The algorithm assumes that 
the X-ray attenuation coefficient of glue is higher than the 
mean wood attenuation coefficient. This is typically the case 
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Fig. 3 (iCT tomoslices, 
a Non-glued material, full-angle 
tomography, b Non-glued, 
limited-angle, c Glued, 
full-angle, d Glued, 
limited-angle 

Abb. 3 pCT Schnittbilder. 
a Nicht verklebtes Holz, 
Vollwinkel Tomographie. 
b Nicht verklebt, 
eingeschrankter Winkelbereich. 
c Verklebt, Vollwinkel. 
d Verklebt, eingeschrankter 
Winkelbereich 
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has highest effective contrast. An optimum performance was 
achieved for A = 2 by substituting \fi a — i± w |A = |/x 
with /Xg and /x fl the mean X-ray attenuation coefficients ob¬ 
served for adhesive and air, respectively. An absolute mini¬ 
mum, with sufficient contrast with respect to and higher 
than the absolute maximum computed for the same position, 
was identified directly as an air gap h a , with jl = fi ext (h a ). 
Other possible outcomes of Eq. 4 correspond to either glued 
material or non-glued material with an air gap of thick¬ 
ness below the resolution of the measurement; in these cases 
jl = ii ext (h g ) was preliminarily stored. 

The final segmentation step transformed the greyscale 

y\ 

image jl(x,z) into a binary image H(x,z) with possi¬ 
ble outcomes h a and h g . For pixels jl = /JL ext (h a ) directly 

A 

H = h a . For pixels jl = i± ext (h g ), a decision on either glued 
or non-glued timber was established using an optimum prob¬ 
abilistic grey level threshold t 0 \ values of jl over t Q were 

A A 

assigned to H = h g and the rest were considered H = h a . 

A single-threshold maximum a posteriori algorithm 
(STMAP) was designed to find t 0 based on training infor¬ 
mation acquired from the samples. This algorithm is a mod¬ 
ification of the classical Bayesian hypothesis test (Van Trees 
2001). The resulting optimum threshold is well defined and 


unique against incomplete estimates of the grey level proba¬ 
bility density functions for glued and non-glued timber. The 
3D joint probability density function of the measured X-ray 

yv 

attenuation coefficient M occurred hypothesis H = {h a , h g ) 

yv 

and estimated hypothesis H = {h a , h g ) as a function of the 
threshold t follows the following expression: 
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PH(h)f^, H (ii\h)U(t -fi) H = h 
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The upper case letters //, H and M denote random vari¬ 
ables and the lower case letters h , jl and jl realisations; 

is the probability density function of the X-ray atten- 


fu 


J M\H 

uation coefficient conditioned on the occurred hypothesis, 


PhQi) are hypotheses probabilities and U(x) is the unit- 
step Heaviside function. The goodness of decision D is a 
binary random variable with outcome “ok” if the estimated 
hypothesis equals the occurred hypothesis and “false” other¬ 
wise. Operating Eq. 5, the probability of a correct decision 
conditioned on the measured X-ray attenuation coefficient 
p D |^(“ok”l/x) is calculated as 


4^ Springer 


























































































610 


Eur. J. Wood Prod. (2011) 69:605-617 


p D \m(“ °k”l£) 


information, leading to the false interpretation that there is 
glue between timber lamellas h g . In the case of Fig. 3c the 
profile C-C has been drawn across EW; a maximum here 
indicates the presence of adhesive in the glue line. 

The limited-angle reconstruction effectively filters out 
the year ring structure and other anatomical wood features 
from the tomograms, resulting in a homogeneous grey level 
in the wood regions (Figs. 3b and d). Comparing A-A with 
its limited-angle version B-B, one observes that LW fea¬ 
tures (2) and (4) are not present in the latter; however, 
the information of the glue line is preserved. A clear grey 
level minimum in B-B is now correctly interpreted as non- 
glued wood h a . In the case of profile C-C, the improvement 
achieved in D-D is not so evident. From several similar ob¬ 
servations it can be concluded that with limited-angle re¬ 
construction, the grey level contrast between adhesive and 
wood and the grey level variability for wood are of the same 
order of magnitude as the values measured in EW regions 
with full-angle reconstruction. The grey level variations in 
the limited-angle case follow a Gaussian distribution. 

Figure 3b shows that false maxima are still possible in 
the limited-angle reconstruction for very thin air gaps be¬ 
tween timber lamellas (6). This effect justifies the necessity 
for the last segmentation step of the algorithm to differen¬ 
tiate between glued regions and thin air gaps. A possible 
explanation is linked to an edge distortion effect; for ex¬ 
ample, in profiles B-B and D-D positive and negative over¬ 
shot is clearly observed in the transitions air/wood, which 
can lead to false maxima detection when the separation be¬ 
tween timber lamellas is below the effective resolution of 
the tomogram. The overshot can be qualitatively explained 
as Gibbs ringing artefact (Archibald and Gelb 2002) intro¬ 
duced by the discontinuity in the Fourier domain between 
reconstructed and non-reconstructed spectral components. 


z(ha,h a ,fr)pH(h a ) + ffi 
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The optimum threshold t Q that maximises the probability of 
correct decision is 

argmax ? po(“ ok”) 


argmax, p D .^(“ak”\p,)dp, 
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arg max, [p H (h g )F^ |ff (t\h g ) 

Ph (.h-a) 

(t\h) = f!_ 0 O fM {H (£\h)djl 


(J) 


F. 


1 M\H 

where F^ H {jX\h g ) and F^ H {jX\h a ) are cumulative prob¬ 
ability distributions of the X-ray attenuation coefficient 
conditioned on the occurred hypotheses. F^^ H (jl\h g ) and 

Fm\h (Al^fl) are estimated by integrating a histogram of 
training pixels for glued and non-glued regions, respectively. 
The main advantage of this procedure is that even if scarce 
or no histogram data is available for certain grey level values 
/x, the cumulative distributions F^ H (jl\h) are well defined 
monotonically increasing functions. Moreover, once t Q is 
available, the implementation of the segmentation process is 
portable to very simple hardware architectures. 


3 Results and discussion 


Figure 3 shows typical tomoslices for well-glued and non- 
glued timber, reconstructed with full-angle (1800 projec¬ 
tions with 0.2° steps) and limited-angle (19 projections, an¬ 
gular range of ±1.8°) jaCT. The full-angle results clearly re¬ 
veal the year ring structure with significantly stronger X-ray 
attenuation in the high density latewood (LW) zone (2) 
than in the lower density early wood (EW) regions. Several 
anatomical features such as axial resin ducts (1), knots (3), 
wood rays (5), and resin pockets (7) can be observed. The 
density was computed from a linear calibration with the 
grey level values /x, using the nominal densities of air 
(1.2 kgm -3 ) and polyurethane adhesive (1100 kgm -3 ) as 
reference values. The resulting wood density values are 

260-340 kgm -3 for EW and 800-1100 kgm -3 for LW. In 

LW areas, /x values of the same magnitude order as the ones 
for adhesive can be observed (2). This can lead to failure 
in the glue line segmentation when the year ring structure 
intersects the glue line (4). Due to the limited spatial res¬ 
olution of the tomograms, the LW information is smeared 
into the glue joint. This effect can be clearly seen in the A-A 
profile, where the EW/LW transition (4) hides the air gap 


3.1 Optimum segmentation of glued and non-glued areas 


Figure 4 illustrates glue line assessment results for a typical 
glued timber sample. The central image is an average optic 
image of the inner lamellas of the open sample. Figures 4a, 
c, e and g are greyscale images /x(v, z) obtained as results 
from the glue line readout algorithm. Generally three differ¬ 
entiated regions can be distinguished in the images, a dark 
grey region R a corresponding to values of pL e xt(h a ) where 
the air gap thickness was large enough to assess lack of glue 
without STMAP segmentation; a confusion region R ? for 
very thin air gaps, where ii e xt(h g ) is plotted; and a well- 
glued region R g where fi ex t(h g ) is also plotted and the X- 
ray attenuation coefficient is largest. Figures 4b, d, f and h 

A 

are binary images H(x,z ) (white: glued, black: non-glued) 
obtained after STMAP segmentation. For all output images 
most glue line pixels in R ? and R g were assigned correctly 
to non-glued and glued material, respectively. 
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Fig. 4 Glue line assessment 
results, a, c, e and g are 
greyscale images jl(x,z) 
calculated according to Table 1 . 
b, d, f and h are the 
corresponding binary images 

H(.x,z) 

Abb. 4 Ergebnisse des 
Klebfugen-Auslesealgorithmus. 
a, c, e und g sind 
Graustufenbilder jl(x,z) gemaB 
Tabelle 1. b, d, f und h sind 
entsprechende 

A 

Binardarstellungen H(x, z) 


a) 


h) 


9 ) 


The reconstruction parameters had a direct impact on 
the final binary images. Figures 4a and b show results for 
full-angle reconstruction (1800 projections with 0.2° steps), 
which resulted in the largest R ? region and worst detec¬ 
tion performance, with some false hits in both R? and R 
The year ring pattern introduces alternating bright (EW) 
and dark (LW) horizontal stripes. This is a consequence 
of a mixture of the X-ray attenuation coefficients of wood 
and glue line for gap thicknesses below the measurement 
resolution. Figures 4c and d show limited-angle results 
(19 projections, angular range of ±1.8°); a remarkable de¬ 
tection improvement is observed, with almost a total ab¬ 
sence of Ri. The transition profile between glued and non- 
glued regions in Fig. 4d is smoother than in Fig. 4b. How¬ 
ever, the streaky features in the transition profile of Fig. 4b 
could not be matched reliably to those measured with op¬ 
tic means. Therefore, it is not clear whether they are re¬ 
ally present in the sample or just artefacts. The smallest de¬ 
tails that were successfully identified in both jaCT and op¬ 
tic images, such as the small adhesive droplet (1) or fea¬ 
tures (2) and (3), have diameters of 5 mm or larger and are 
detected with both full-angle and limited-angle reconstruc¬ 
tion. 


CNR = 20 log 


( 8 ) 


.V 


/x 


where /n g is the effective X-ray attenuation coefficient of the 
adhesive detected in the glue line; /Ji e g equals jl in R g . The 
CNR decreases for large grey level variations in wood or 
for small fi g , with the latter associated to glue line thick¬ 
nesses under the resolution limits of the measurement. Fig¬ 
ure 5 shows CNR statistics as a function of the angular range 
reconstructed. The CNR decreases monotonically with in¬ 
creasing number of reconstructed projections, for an angu¬ 
lar range of about ±2° the observed CNR is between 20 and 
25 dB, for full-angle reconstruction it varies between 6 and 
20 dB (smaller CNR values than 6 dB are assigned to R? 
according to k in Eq. 4). The CNR shows a sharp increase 
between 1701 (angular range of ±170°) and 1800 (±180°) 
projections. In the latter case the glue line plane is irradiated 
from both sides of the sample, thus causing an averaging by 
a factor of 2 of the glue line information that theoretically 
reduces the random noise level by 3 dB. Best CNR is ob¬ 
tained for radiography (single projection across the glue line 
plane), which expresses the inherent trade-off in the recon¬ 
struction between contrast (CNR) and spatial resolution in 
the bonding plane. In general, 9 projections (angular range 
of ±0.8°) were the minimum required to detect glue line 
features of one tenth of the sample thickness, best resolution 
was obtained for 19 projections (±1.8°). Significant noise 


8 ' 


A quantitative assessment of the adhesive detectability 
with respect to the heterogeneity of the surrounding wood 
material was performed in terms of a logarithmic contrast- 
to-noise ratio (CNR) defined as follows: 
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Fig. 5 Contrast to Noise Ratio 
(CNR) as function of limited 
angular range. Dashed lines are 
mean values for each glued 
timber sample measured, for one 
of them {solid line ) confidence 
intervals are also plotted 
Abb. 5 CNR in Abhangigkeit 

des eingeschrankten 
Winkelbereichs. Gestrichelte 
Linien sind Mittelwerte fur jede 
gemessene verklebte Holzprobe; 
Konhdenzintervalle sind fur 
eine davon (ausgezogene Linie ) 
auch gezeichnet 
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Table 1 pCT parameters and computational costs corresponding to the results of Fig. 4 
Tab. 1 (iCT BestimmungsgroBen und Rechenaufwand der Ergebnisse in Abb. 4 


Fig. 4a-b 


Fig. 4c-d 


Fig. 4e-f 


Fig. 4g-h 


Angle step (°) 

No. projections 

Frame pixel size Ui x vj (pm 2 ) 
Reconstructed voxel v; x yj x zk (prn 3 ) 
CNR (dB) 

FOV e (mm) 

Memory space X-ray frames 3 
Memory space 3D tomogram 3 
X-ray data acquisition time 
Backprojection time 


0.2 


0.2 


0.2 


0.4 


1800 


19 


19 


9 


41 x 41 

164 x 41 x 164 


41 x 41 

164 x 41 x 164 


41 x 820 
820 x 41 x 820 


82 x 1640 
1640 x 82 x 1640 


16 


21 


25 


25 


112 


18 


18 


17 


14.4 GB 
180 MB 

120 min 
760 min 


25 MB 
180 MB 


1.3 MB 

7.3 MB 


150 KB 
0.92 MB 


1.3 min 
7.9 min 


1.3 min 


36 s 


15 s 


1.2 s 


Data stored with 16 bits per sample. A volume of 100 x 10 x 100 mm 3 was reconstructed 


pixel size is 41 x 820 (dm 2 ). A clear noise reduction is ob¬ 
served in comparison to Fig. 4d, with a CNR improvement 
of 4 dB. However, it is also observed that the number of 
pixels belonging to R? increases, this effect can be inter¬ 
preted as a resolution loss in y which results in a smearing 
of the X-ray attenuation coefficient in the glue line. Bin¬ 
ning in u had a direct impact on the glue line readout; for 
binning factors larger than 2 the assessment was no longer 
possible. The glue line assessment performance was not af¬ 
fected by cropping the detector frame in u to FOV e cal¬ 
culated with Eq. 2. Finally, the combination of several in¬ 
formation reduction techniques was investigated. Figures 4g 
and h were calculated with 9 projections spaced 0.4°, a bin¬ 
ning factor of 2 in u and a binning factor of 40 in v (sample- 
related pixel size of 82 x 1640 (dm 2 ). Even in this case, the 
adhesive droplet (1) and feature (3) were clearly recognis¬ 
able. 


degradation with no further resolution improvement was ob¬ 
served above 51 projections (±5°). 


3.2 Minimisation of reconstruction information 


The reconstruction information was minimised by enlarging 
the angular steps between projections P$> and by introduc¬ 
ing an anisotropic binning factor in both frames and recon¬ 
structed tomograms. Angular steps of 0.15° did not signif¬ 
icantly improve the glue line assessment in comparison to 
steps of 0.2°. With 0.4° steps the minimum number of pro¬ 
jections could be reduced from 9 to 7 with a 3 dB CNR 
decrease; steps of 0.8° or larger degraded the binary images. 
A reliable assessment was not possible using asymmetric an¬ 
gular projections. Binning in v decreased the resolution in z 
but also improved the CNR due to the averaging applied to 
the X-ray detector pixels. Figures 4e and f show the out¬ 
come of a reconstruction with 19 projections (angular range 
±1.8°) binning with a factor of 20 in v (the sample-related 


Table 1 summarises the reconstruction settings used 
to compute Fig. 4 together with time and memory costs. 
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The backprojection has a computational complexity of 
0(N x N y N z N&), where Ni are the number of pixels in each 
dimension. The acquisition time was proportional to N$> and 
the exposition time of the X-ray detector. The total memory 
space for a 100 x 100 x 10 mm 3 volume was reduced by a 
factor of 14000 between full-angle (Fig. 4a-b) and limited- 
angle reconstruction with multiple information reductions 
(Fig. 4g-h). The acquisition time was reduced by 200 and 
the backprojection time by 38000. The glue line assessment 
of Fig. 4g-h was implemented in less than 40 s. It is pos¬ 
sible to significantly reduce the computation times shown 
in Table 1 ; for example, by using a low level programming 
language or by parallelising the discussed algorithms. New 
developments in Graphic Processing Units (GPUs) allow for 
low-cost hardware parallelisation of the FDK algorithm with 
a time reduction by a factor of 80 with respect to a single- 
threaded CPU (Sharp et al. 2007). Combining these optimi¬ 
sations, a real time implementation of the glue line assess¬ 
ment is feasible. The data acquisition times could be reduced 
by at least one order of magnitude by using a directional tar¬ 
get X-ray source (for example XS160D, GE Sensing & In¬ 
spection Technologies GmbH, Wunstorf, Germany), which 
allows for larger tube currents, therefore minimising the ex¬ 
position times of the detector. 


of this setup was chosen in order to achieve sufficient field 
of view for the full-angle measurements. 

Figure 6 shows the stereomicroscopy measurements of 
the glue line thickness in the vicinity of the glued/non-glued 
transition and along a typical adhesively bonded glue line. 
The gap thickness readout was implemented with a similar 
algorithm to the one described in the next section. A total of 
14 transitions were analysed. The results are summarised in 
Table 2. The mean air gap thickness in an interval of 1 mm 
from the transition was 15 (dm, the adhesively bonded glue 
line thickness was 50 (am. Typical adhesive joint thicknesses 
in a 10 mm interval separated by at least 10 mm from the 
transition oscillated between 100 and 200 jam. The air gap 
was only thicker than 150 jam at distances from the transi¬ 
tion that were greater than 20 mm. The position error be- 
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3.3 Spatial resolution and detectability limits 
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urn 


The effective full-angle spatial resolution across the glue 
line was measured with a self-developed algorithm that cal¬ 
culated the distance between 10% and 90% level cuts in the 
air/wood transitions (Fig. 3c). Statistics were calculated for 
the full air/wood transition area of 10 lamellas, with an ex¬ 
cellent fit to a Gaussian distribution with a mean of 140 (dm 
and a std. deviation of 40 (dm. The main resolution limit is 
given by the effective pixel size of the X-ray detector (ac¬ 
cording to the datasheet 105 (dm sample-related pixel size). 
The unsharpness of the X-ray beam is less than 20 |dm due to 
the large distance between source and sample. The effective 
sample-related pixel size can be reduced in about 25% by 
placing the sample closer to the X-ray tube. The geometry 
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Fig. 6 Glue line thickness characterisation using stereomicroscopy, 
a Transition between glued and non-glued timber, b Typical adhesive 
joint thickness profile 

Abb. 6 Charakterisierung der Dicke der Klebfuge mittels Stereomi- 

• • 

kroskopie. a Ubergang zwischen verklebtem und nicht verklebtem 
Holz. b Dickenprohl einer typischen Leimschicht 


Table 2 pCT spatial resolution 
across the bonding plane in 
comparison to stereomicroscopy 
results 

Tab. 2 Raumliche Auflosung 
der pCT quer durch die 
Leimschicht im Vergleich zu 
Ergebnissen mittels 
Stereomikroskopie 


Mean estimated value (pm) 


Parameter 


pCT sample-related pixel size 
pCT effective sample-related pixel size 
Stereomicroscopy effective pixel size 
Glue line thickness (1 mm from transition) 

Air gap thickness (1 mm from transition) 

Glue line thickness (> 10 mm from transition) 

Air gap thickness (>20 mm from transition) 

Transition positioning error between stereomicroscopy and pCT 


41 


105 


35-60 

10-20 

100-200 


150 


1000-2000 
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Fig. 7 (iCT glue line topology 
assessment, a Gap thickness 
map. b Open board photograph, 
c Tomoslice for profile A-A. 
d Comparison between pCT and 
optic air gap thickness 
determination 
Abb. 7 Messung der 
Klebfugendicke mittels pCT. 
a Abbildung der Spaltdicke. 
b Photographic der 
aufgebrochenen Lamellen. 
c Schnittbild fur Prohl A-A. 
d Vergleich zwischen pCT und 
optisch ermittelter 
Luftspaltdicke 
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tween transitions measured with jaCT and stereomicroscopy 
was typically smaller than 2 mm. From these results it can 
be deduced that the |dCT method detects the presence or 
absence of adhesive in glue lines of thicknesses down to 
half the effective sample-related pixel size. In the vicinity 
of the transition, the thin air gap between timber lamellas 
falls well below the resolution of the tomograms and cannot 
be detected. However, the presence or absence of adhesive 
in these confusion regions Ri can still be assessed with the 
segmentation algorithms previously described. 


Figure 7a shows a typical jaCT gap thickness topology 
assessment, an optic image of the open sample is shown in 
Fig. 7b. An angular range of ±2.6° (35 projections) was cho¬ 
sen. A small knot on the board surface (1) introduced large 
X-ray attenuation coefficients in the jaCT images, which 
could be misinterpreted as glued material. However, the gap 
topology assessment reveals the surface topology of the de¬ 
fect, providing information about its nature. The air gap 
thickness (2) and the adhesive joint (3) topology could also 
be investigated. Due to the presence of an aluminium spacer 
at middle length of the sample (Fig. 1) it was observed that 
the wood lamellas were slightly curved, being in close con¬ 
tact at the centre of the width and showing an opening (2) 
at the edges (Fig. 7c). Air gap thickness profiles measured 
optically at the edges of the sample and gap thickness pro¬ 
files evaluated from jaCT tomograms (Fig. 7d) were com¬ 
pared. Two optic measurements were performed immedi¬ 
ately before and after jaCT measurement, a 100 jam variation 
was observed between both of them. This change was most 
probably due to the manipulation of the sample and creep 
deformation due to the tension induced by the aluminium 
spacer. The jaCT results were in close agreement with the 
optic data for air gap thicknesses larger than 150 jam, which 
corresponds to the effective resolution limit of the jaCT mea¬ 
surement. 


3.4 Gap topology assessment 


Next, the extension of the glue line readout algorithm for gap 
topology assessment between wood lamellas is discussed. 
The digital resolution of the tomograms across the glue line 
(y direction) was increased by a Fourier domain interpola¬ 
tion. The gap thickness between glued timber lamellas was 
computed for each position (a, z) of the bonding plane by 
recording the width of the positive or negative peak observed 
in the glue line, thresholded with a cut value /x CM? (a, z,h): 


t 


8 


Hcutihg) — fi w + 


[f^ext (hg) Mw] 


100 


( 9 ) 
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H'cut (Ha) — 


f^ext (ha)] 


[m 


XJO 


XJ0 


100 


where t g and t a are percentage thresholds for glued and non- 
glued material, respectively. They depend on the effective 
contrast of adhesive and air with respect to wood, related 

\h < | l^ext(hg) 

^cut(hg )I itfol- 

- G/100)/(l — /g/100). Best results were 
50% and t 


4 Results on glulam lamellas of thicknesses used in 
construction 


to X with Eq. 4. Rewriting \ii e xt(h a ) 

l^w \ &S \f^ext (ha) ^cut (ha) I < I f^ext(hg) 

lows that X = (1 
observed for t 


Finally, glue line assessment results for glued timber lamel¬ 
las of the thicknesses typically found in glulam construc¬ 
tions (EN 386:2001) are shown. Two lamellas of 250 x 
170 x 32.5 mm 3 were glued together to form a glulam beam 


75%. 
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Fig. 8 (iCT assessment of 
glulam. a Full-angle (top left ) 
and limited-angle (top right ) 
tomoslices, b and c are 
corresponding jl(x,z). 

A 

d H(x,z ) computed from (b). 
e Air-coupled ultrasound results 
Abb. 8 pCT Messungen an 
Brettschichtholz. a Schnittbilder 
fur Vollwinkel (oben links ) und 
eingeschrankten Winkelbereich 
{oben rechts ). b und c sind 
entsprechende jl(x, z). 

A 

d H(x,z ) berechnet aus (b). 
e Ergebnisse mit 
luftgekoppeltem Ultraschall 
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of 250 x 170 x 65 mm 3 . Adhesive was applied to the en¬ 
tire sample except for an 80 x 80 mm 2 region. A tube volt¬ 
age of 80 kV with a nominal power of 15 W was used, 
the detector pixel was hardware binned to 100 jam for a 
sample-related pixel size of 87 jam and an exposition time 
of 0.5 s was chosen. 1800 projection frames were recorded 
with 0.2° steps. Figure 8a shows reconstructed tomoslices 
using 1800 and 19 projections (no frame binning, voxel size 
of 697 x 87 x 697 jam 3 ). The region (1) fell outside the field- 
of-view only in the full-angle case; (5) was reconstructed 
outside the glued timber sample. The year rings were essen¬ 
tially parallel to the adhesive joint, which introduced strong 
greyscale variations across the glue line. In contrast to Fig. 3, 
the year ring pattern of both lamellas did not match in most 
of the bonding plane. In the full-angle case, the EW/LW 
transitions present in the glue joint hided the adhesive in¬ 
formation (profile A-A) and the glue line assessment failed 
to detect the lamination faults (Fig. 8b). The limited-angle 
reconstruction revealed the presence of adhesive (2) in the 
glue line (profile B-B) and filtered out to a great extent the 
year ring structure except for EW/LW transitions (3) highly 
parallel to the glue line (angle smaller than 1.8°). The re¬ 
maining wood structure did not compromise the imaging of 
the bonding plane if a small assessment interval was used; 

2 mm was sufficient. Figures 8c and d show limited- 
angle glue line assessment results. The non-glued region 
was successfully identified and the EW/LW transitions in¬ 
troduced only a marginal number of noisy pixels (4) in the 
binary image. However, the measured non-glued region was 
smaller (about 30 x 70 mm 2 ) than the defect area defined 


during the manufacturing of the sample (dashed square). 
A probable reason is that the adhesive may have penetrated 
into the defect area during the pressing of the glued tim¬ 
ber lamellas. The jaCT images were compared to an air- 
coupled ultrasound (ACU) measurement of the same sample 
(Fig. 8e) according to Sanabria et al. (2010b). A decrease 
of ultrasonic amplitude revealed a reduced defect area that 
closely matched the position and extension observed with 
the jaCT assessment. 


5 Conclusion 


It was shown that Limited-angle Microfocus X-ray Com¬ 
puted Tomography is an interesting alternative for glue line 
assessment of timber constructions; inspection of pairs of 
100 and 170 mm wide glued timber lamellas is possible. 
A simple limited-angle FDK reconstruction with 7 to 19 
projection frames in an angular range of ±0.6° to ±1.8° 
was sufficient to segment the bonding plane into glued and 
non-glued areas with a spatial resolution of 5 mm. The 
small angular ranges allowed for fields of view of less than 
20 mm, making the measurement independent of the height 
and length of the inspected beam. Moreover, the limited- 
angle reconstruction efficiently filters out undesired wood 
structure while preserving the glue line information. Accu¬ 
rate segmentation of the glue line was possible in all in¬ 
spected samples, even if similar X-ray attenuation was ob¬ 
served for adhesive and latewood regions. The measurement 
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time was reduced to 40 s for a 100 x 100 mm 2 bonding plane 
by reconstructing only a reduced number of anisotropically 
binned projections. A detailed analysis of the resolution of 
the (dCT tomograms together with stereomicroscopy mea¬ 
surements of the glue line thickness revealed that, in the 
proximity of the transitions between glued and non-glued ar¬ 
eas, the air gap separation between timber lamellas fell well 
below the 105 (dm resolution of the (dCT measurement. Edge 
distortion effects in the air/wood transition can lead to a con¬ 
fusion region in which glued material is falsely detected. 
A STMAP algorithm was developed and allowed for a ro¬ 
bust probabilistic discrimination of glued material from thin 
air gaps. The complete measurement procedure was able to 
detect the presence or the lack of adhesive down to glue 
lines of 50 (dm thickness. Delaminations and air gap thick¬ 
ness topology could be assessed for separations larger than 
150 (dm between timber lamellas. These limits did not com¬ 
promise the detectability of lamination faults for thin glue 
lines of 100 to 200 (dm obtained with hydraulic pressing. 

Future research work should extend the implementation 
to an industrial environment. The described method can be 
applied separately to individual laminations; therefore, its 
extension to multilayered glulam assessment is immediate. 
The geometric arrangement described should be improved 
to allow for practical positioning of large samples with re¬ 
spect to the X-ray tube and detector. A real time implemen¬ 
tation of the measurement procedure seems feasible, with 
the exposition times in the detector being reduced by using 
a directional target source. New developments in Graphic 
Processing Units (GPUs) should allow for low-cost real time 
implementation of the glue line assessment algorithms. 


Dill-Langer G, Aicher S, Bemauer W (2005) Reflection measurements 

at timber glue-lines by means of ultrasound shear waves. Otto- 

Graf-J 16:273-283 

Dunky M, Niemz P (2002) Holzwerkstoffe und Leime: Technologie 

und Einflussfaktoren. Springer, Berlin 
EN 386:2001 Glued laminated timber—Performance requirements and 

minimum production requirements 

EN 391:2001 Glued laminated timber—Delamination test of glue lines 
EN 14080:2005 Timber structures—Glued laminated timber— 

Requirements 

Entacher K, Planitzer D, Uhl A (2007) Towards an automated genera¬ 
tion of tree ring profiles from CT-images. In: Proc of the 5th inter¬ 
national symposium on image and signal processing and analysis, 
New York, NY, USA, pp 174-179 

Espinoza GR, Hernandez R, Condal A, Verret D, Beauregard R (2005) 

Exploration of the physical properties of internal characteristics of 
sugar maple logs and relationships with CT images. Wood Fiber 

Sci 37(4):591-604 

Feldkamp LA, Davis LC, Kress JW (1984) Practical cone-beam algo¬ 
rithm. J Opt Soc Am A, Opt Image Sci Vis 1 (6):612—619 
Hasenstab A (2006) Integritatspriifung von Holz mit dem zerstorungs- 

freien Ultraschallechoverfahren. Dissertation, Technische Univer- 
sitat Berlin 

Hu LJ, Gagnon S (2007) X-ray-based scanning technique for non¬ 
destructive evaluation of finger-joint strength. In: Proc of the 15th 
international symposium on NDT of wood, Deluth, MN, USA 
Kabir MF, Schmoldt DL, Schafer ME (2002) Time domain ultrasonic 

signal characterization for defects in thin unsurfaced hardwood 

lumber. Wood Fiber Sci 34(1): 165-182 
Kak AC, Slaney M (1988) Principles of computerized tomographic 

imaging. IEEE Press, New York 

Longuetaud F, Mothe F, Leban JM (2007) Automatic detection of the 

heartwood/sapwood boundary within Norway spruce (Picea abies 
(L) Karst) logs by means of CT images. Comput Electron Agric 

58:100-111 

Maeva E, Severina I, Bondarenko S, Chapman G, O’Neill B, Severin 

F, Maev RG (2004) Acoustical methods for the investigation of 
adhesively bonded structures: a review. Can J Phys 82(12):981— 
1025 

Mannes D, Sonderegger W, Hering S, Lehmann E, Niemz P (2009a) 

Non-destructive determination and quantification of diffusion 
processes in wood by means of neutron imaging. Holzforschung 

63(5):589-596 

Mannes D, Niemz P, Lehmann E (2009b) Tomographic investiga¬ 
tions of wood from macroscopic to microscopic scale. Wood Res 

54(3):33-44 

Morigi MP, Casali F, Bettuzzi M, Bianconi D, Brancaccio R, Comac- 

chia S, Pasini A, Rossi A, Aldrovandi A, Cauzzi D (2007) CT 
investigation of two paintings on wood tables by Gentile da Fab¬ 
riano. Nucl Instrum Methods Phys Res, Sect A, Accel Spectrom 
Detect Assoc Equip 580(l):735-738 
Niemz P, Kucera LJ, Flisch A, Blaser E (1997) Application of comput- 

ertomography (CT) on wood defects and decay. Holz Roh- Werkst 

55:279-280 

Niemz P, Poblete H, Cuevas H, Flisch A (1999) Klebfugenfestigkeit 

von keilverzinktem. Holzforsch Holzverwert 51:79-81 
Oja J (2000) Evaluation of knot parameters measured automatically in 

CT-images of Norway spruce (Picea abies (L) Karst). Holz Roh- 

Werkst 58(5):375-379 

Osterloh K, Raedel C, Zscherpel U, Meinel D, Ewert U, Buecherl T, 

Hasenstab A (2008) Fast neutron radiography and tomography of 

wood. Insight 50(6):307-311 

Rinnhofer A, Petutschnigg A, Andreu JP (2003) Internal log scanning 

for optimizing breakdown. Comput Electron Agric 41 (1-3):7-21 
Sanabria SJ, Mueller C, Neuenschwander J, Niemz P, Sennhauser U 

(2010a) Air-coupled ultrasound as an accurate and reproducible 


Acknowledgement This research has been supported by the Swiss 
National Science Foundation under contract 200021-115920. The au¬ 
thors acknowledge the work of Fabian Binkert, Danniel Gomes and 
Dominik Machler for the optical characterisation of the samples and 
analysis of pCT data, and the tutorials of Dr. Iwan Jerjen and Dr. Rene 
Mooser in CT reconstruction software. 


References 


Archibald R, Gelb A (2002) A method to reduce the Gibbs ringing ar¬ 
tifact in MRI scans while keeping tissue boundary integrity. IEEE 
Trans Med Imaging 21(4): 305-319 
Berglind H, Dillenz A (2003) Detection of glue deficiency in laminated 

wood with pulse thermography. J Wood Sci 49(3):216-220 
Bhandarkar SM, Luo XZ, Daniels R, Tollner EW (2006) A novel 

feature-based tracking approach to the detection, localization, and 
3-D reconstruction of internal defects in hardwood logs using 
computer tomography. Pattern Anal Appl 9(2-3): 155-175 
Bucur V (2003) Nondestructive characterization and imaging of wood. 

Springer, Berlin 

Choi MY, Park JH, Kim WT, Kang KS (2008) Detection of delami¬ 
nation defect inside timber by sonic IR. In: Proc of the Society 
of Photo-optical Instrumentation Engineers (SPIE), Orlando, FL, 

USA, 6939 U306-U309 


4^1 Springer 







Eur. J. Wood Prod. (2011) 69:605-617 


617 


method for bonding assessment of glued timber. Wood Sci Tech- 
nol. doi:10.1007/s00226-010-0357-z 
Sanabria SJ, Neuenschwander J, Niemz P, Sennhauser U (2010b) 

Structural health monitoring of glued laminated timber with a 
novel air-coupled ultrasound method. In: Proc of the 11th World 
conference on timber engineering, Riva del Garda, Trentino, Italy 
Scheepers G, Moren T, Rypstra T (2007) Liquid water flow in Pinus 

radiata during drying. Holz Roh- Werkst 65(4):275-283 
Schmoldt DL, He J, Abbott AL (2000) Automated labeling of log fea¬ 
tures in CT imagery of multiple hardwood species. Wood Fiber 

Sci 32(3):287-300 

Sepulveda P, Oja J, Gronlund A (2002) Predicting spiral grain by com¬ 
puted tomography of Norway spruce. J Wood Sci 48(6):479-483 


Sharp GC, Kandasamy N, Singh H, Folkert M (2007) GPU-based 

streaming architectures for fast cone-beam CT image recon¬ 
struction and demons deformable registration. Phys Med Biol 

52(19):5771-5783 

Sirr SA, Waddle JR (1999) Use of CT in detection of internal dam¬ 
age and repair and determination of authenticity in high-quality 
bowed stringed instruments. Radiographics 19(3):639-646 
Solodov I, Pfleiderer K, Busse G (2004) Nondestructive characteriza¬ 
tion of wood by monitoring of local elastic anisotropy and dy¬ 
namic nonlinearity. Holzforschung 58(5):504-510 
Van Trees HJ (2001) Detection, estimation and modulation theory, 

Part I. Wiley, New York 


4^) Springer 



